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Density  functional  theory  is  applied  to  the  computation  of  the  adsorption  energy  (AEads)  for  a  series  of 
molecules  on  y-A^Os.  Three  different  cluster  models  are  used  to  represent  the  y-A^Oa  surface.  The  mol¬ 
ecules  of  interest  all  contain  a  phosphonyl  (P=0)  functional  group  and  adsorb  via  formation  of  a  donor 
bond  between  the  0  atom  and  a  threefold-coordinated  tetrahedral  A1  [Al(Td)]  surface  site.  The  highest 
occupied  molecular  orbital  (HOMO)  of  the  free  molecule  is,  in  all  cases,  composed  largely  of  non-bonding 
orbitals  on  the  0  atom  of  the  P=0  group.  The  empty  “dangling  orbital”  on  the  coordinatively-unsaturated 
Al(Tci)  site  constitutes  a  surface  state.  A  linear  relationship  is  found  between  Afads  and  the  difference 
between  the  orbital  energies  of  the  molecular  HOMO  (sh)  and  the  surface  state  (es)-  Trends  in  AEads  for 
different  molecules  can  then  be  understood  in  terms  of  variations  in  6h-  Likewise,  differences  in  Afads 
for  various  cluster  models  can  be  explained  by  considering  the  differences  in  the  predicted  £s  values.  A 
further  evaluation  of  the  cluster  models  is  presented  by  comparing  results  for  the  physisorption  of 
H2O  or  CO  with  those  obtained  from  two-dimensionally-periodic  slab  models.  When  differences  in  en 
and  £s  4re  accounted  for,  the  various  models  and  computational  procedures  are  seen  to  yield  essentially 
equivalent  results  for  adsorption  of  the  molecules  considered.  These  results  are  thought  to  constitute  a 
useful  conceptual  tool  for  rationalizing  AEads  values  for  different  molecules  and  cluster  models. 
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1.  Introduction 

Many  highly  toxic  chemicals  involve  the  phosphonyl  (P=0) 
functional  group.  These  include  herbicides  and  pesticides  as  well 
as  chemical  warfare  agents  (CWAs)  such  as  Tabun  (GA),  Sarin 
(GB),  Soman  (GD)  and  VX.  The  adsorption  of  these  species  on  var¬ 
ious  substrates  is  an  important  issue  in  “agent  fate”,  i.e.,  the  stabil¬ 
ity  of  the  adsorbate  with  respect  to  desorption,  decomposition  and 
reaction  with  environmental  species.  In  both  experiment  [1,2]  and 
theory  [3-6],  adsorption  on  the  surface  of  a  Lewis-acid  material 
such  as  OH-free  y-A^Os  occurs  via  dative  bond  formation  involving 
the  phosphonyl  0  atom  and  a  coordinatively-unsaturated  cation 
surface  site.  Non-bonding  orbitals  (NBOs)  on  the  former  comprise 
the  highest-occupied  molecular  orbital  (HOMO)  of  the  free  mole¬ 
cule.  Unoccupied  orbitals  localized  mainly  on  the  surface  cations 
form  a  band  of  surface  states  on  the  oxide. 

One  purpose  of  this  work  is  to  explore  the  correlation  between 
the  adsorption  energy  (AEads)  and  the  energy  of  the  HOMO  relative 
to  that  of  the  surface  state.  This  involves  the  “frontier  orbitals”  con¬ 
cept  of  chemical  bonding  [7-9]  which,  to  our  knowledge,  has  not 
previously  been  applied  to  a  quantitative  description  of  adsorption 
in  these  systems.  The  molecules  studied  (Fig.  1 )  were  chosen  to  ex¬ 
hibit  a  fairly  wide  range  of  ionization  potentials  (IPs,  see  below). 
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These  include  trichlorophosphine  oxide  (TCPO),  dimethyl  methyl- 
phosphonate  (DMMP),  Sarin  and  trimethylphosphine  oxide 
(TMPO).  They  were  also  chosen  to  be  fairly  small  so  that  steric  ef¬ 
fects  would  not  seriously  complicate  the  interpretation  of  AEads- 
DMMP  is  also  of  interest  since  it  is  frequently  used  as  a  safe  simu¬ 
lant  for  Sarin.  Another  purpose  is  to  understand  how  the  choice  of 
cluster  model  affects  the  final  results  for  the  adsorption  of  these 
species.  This  study  involves  several  issues  and  computational  ap¬ 
proaches.  The  properties  of  the  free  molecules,  the  surface-state 
structure  of  the  bare  semi-infinite  y-AlOs  surface  and  of  different 
cluster  models  and  the  interaction  of  the  molecules  with  the  clus¬ 
ters  all  have  to  be  considered  in  unison. 


2.  Computational  details 

2.1.  Molecular  properties 

The  free  molecules,  and  also  various  clusters  designed  to  model 
the  reactive  site  on  the  y-A^Os  surface,  were  treated  using  density 
functional  theory  (DFT)  as  implemented  in  the  Gaussian  03  [10] 
and  the  Amsterdam  density  functional  (ADF)  2006.01  [11]  suites  of 
programs.  The  former  calculations  employed  the  B3LYP  hybrid  func¬ 
tional  with  6-31  lG(d)  basis  sets,  “tight”  convergence  criteria  in 
geometry  optimization  and  an  “ultrafine”  integration  grid.  The  ADF 
calculations  used  the  Vosko-Wilk-Nusair  (VWN)  form  of  the  local 
density  approximation  (LDA)  and  the  Perdew-Burke-Ernzerhoff 
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(d)  TMPO 


Fig.  1.  Molecular  structures  for  TCPO,  DMMP,  Sarin  and  TMPO.  The  Sarin  enantio¬ 
mer  was  arbitrarily  chosen. 

(PBE)  form  of  the  generalized  gradient  approximation  (GGA)  ex¬ 
change-correlation  (XC)  potential.  The  ADF  basis  sets,  designated 
“TZP”,  were  Slater-type  functions  of  triple-zeta  quality  with  a  single 
p  polarization  shell  for  H  and  a  single  d  polarization  shell  for  all  other 
atoms.  These  are  all-electron  basis  sets  with  no  “frozen  cores”  (cf. 
Ref.  [11]).  The  criteria  for  integration  accuracy  and  geometry  conver¬ 
gence  (10“®Hartree  total  energy;  1  x  10“^  Hartree/A  gradient; 
1  X  10“^  A  displacement)  were  tighter  than  the  ADF  defaults  (in  or¬ 
der  to  facilitate  the  computation  of  vibrational  frequencies)  but  loos¬ 
er  than  those  used  in  the  Gaussian  calculations. 

Vertical  ionization  potentials  (VIPs)  were  calculated  using  the 
A(SCF)  method  [12].  Here  the  ground-state  geometry  is  optimized 
and  the  total  energy  obtained.  Then  an  electron  is  removed  from 
the  HOMO  and  the  total  energy  computed,  in  a  spin-unrestricted 
calculation,  with  the  nuclei  fixed  in  the  ground-state  positions. 
The  energy  difference  is  the  VIP  which  can  be  compared  with  data 
from  photoelectron  spectroscopy  (PES).  Adiabatic  ionization 
potentials  (AlPs)  were  obtained  similarly  except  that  the  cation 
geometry  was  allowed  to  relax.  The  AlP  difference  in  energies 
was  adjusted  for  the  difference  in  the  total  vibrational  zero-point 
energy  (ZPE)  which  amounts  to  a  correction  of  ^70  meV  for  the 
present  molecules.  Stability  tests  [13]  were  performed  for  the 
B3LYP  calculations  and  the  cation  wavefunction  was  reoptimized 
in  one  case  (TMPO)  where  an  instability  was  found  (see  below). 
For  B3LYP/6-311G(d)  calculations,  the  harmonic  vibrational  fre¬ 
quencies  were  scaled  by  a  factor  of  0.9669  following  results  [14] 
for  B3LYP/6-311G(d,p).  No  correction  was  applied  for  the 
(VWN  +  PBE)/TZP  calculations  since  the  scaling  factor  in  this  case 
is  estimated  [14]  to  be  close  to  unity. 

2.2.  Bulk  and  two-dimensionally-periodic  slab  models  for  y-Al203 

Calculations  for  bulk  y-A^Os  and  for  a  two-dimensionally  peri¬ 
odic  slab  (2-DPS)  were  done  using  the  Crystal  03  suite  of  programs 
[15,16]  in  order  to  study  surface  states.  The  B3LYP  hybrid  func¬ 
tional  was  used  due  to  its  demonstrated  [17]  ability  to  give  band 
gaps  in  good  agreement  with  experiment  for  a  wide  range  of  mate¬ 
rials  including  AI2O3.  All-electron  Gaussian  basis  sets  [15]  for  A1 
and  0,  designated  86-21dlG  (or  86-21G(d))  and  8-411dllG  (or 
8-411G(2d)),  respectively,  were  used  without  modification  to  con¬ 
struct  Bloch  functions.  The  /<-point  sampling  used  an  8x8x8 
(8x8x1)  grid  for  the  bulk  (2-DPS)  and  an  “extra-large” 


(75,974)  DFT  integration  grid.  In  Crystal,  truncation  of  the  sums 
of  Coulomb  and  exchange  terms  in  the  Fock  matrix  is  determined 
by  five  overlap  criteria  (T1-T5)  [15,16].  These  were  set  at  10“^  for 
T1-T4  and  10“^^  forT5.  In  the  multipolar  expansion  zone  [15,16]  a 
maximum  order  of  L  =  6  was  used. 

The  y-Al203  model  used  here  is  the  defective-spinel  structure 
given  by  Pinto  et  al.  [18]  on  the  basis  of  plane-wave  pseudopoten¬ 
tial  LDA  calculations.  The  bulk-  (and  thus  the  surface-)  structure  of 
y-Al203  is  a  subject  of  controversy  [18-29],  the  essential  aspects  of 
which  are  summarized  elsewhere  [30-32].  Realistic  models  of 
high-surface-area  y-Al203  powders,  even  in  the  absence  of  hydrox- 
ylation,  may  involve  a  wide  variety  of  different  Al  and  0  surface 
sites  as  discussed  in,  for  example.  Ref.  [33].  The  choice  of  bulk 
model  and  surface  plane  determines  the  distribution  and  identity 
of  potential  reactive  sites.  For  example,  the  (100)  surface  of  the 
non-spinel  model  [22]  exposes  three  inequivalent  fivefold-coordi¬ 
nated  Al(Oh)  sites,  each  with  a  different  chemical  character.  The 
present  work  focuses  on  one  particular  type  of  surface  site; 
namely,  the  threefold-coordinated  Al(Td).  This  is  most  easily  mod¬ 
eled  (see  below)  using  the  (1 1  l)a  surface  of  the  defective-spinel 
structure  [18]  but  also  occurs  on  the  (110)  surface  of  the  non¬ 
spinel  model  [22].  Hence,  the  present  work  is  primarily  a  study 
of  adsorption  at  threefold-coordinated  Al(Td)  sites  rather  than  on 
a  particular  surface  plane  of  a  specific  bulk-lattice  model.  It  will 
be  shown  below  that  the  essential  difference  among  these  models 
lies  in  the  energies  of  the  empty  surface  states  derived  from 
cations. 

In  the  present  work,  the  lattice  constants  were  kept  fixed  at  the 
optimized  bulk  values  [18]  while  the  atom  positions  were  allowed 
to  vary  in  a  geometry  optimization  of  the  bulk  structure.  Only  very 
small  differences  (<0.002 )  were  found  between  the  fractional  atom 
positions  and  those  of  the  initial  structure  [18].  A  (1  x  1)  slab  unit 
cell  was  used  with  a  thickness  of  31  atomic  layers  (equivalent  to 
one  crystallographic  unit  cell)  and  the  (1 1  l)a  plane  as  the  surface 
layer.  The  slab  unit  cell  was  symmetric  about  the  mid-plane  so  that 
no  dipole  potential  was  generated.  The  (lll)a,  which  has  been 
identified  [18]  as  the  most  stable  surface  of  a  bulk  single  crystal, 
consists  of  singly-unsaturated  tetrahedral  Al(Td)  sites  and  two 
types  of  0  sites.  One  is  threefold-coordinated  in  the  bulk,  as  a  result 
of  being  adjacent  to  a  vacant  octahedral  Al(Oh)  spinel  site,  and  re¬ 
mains  so  at  the  surface.  This  0  has  one  Al(Td)  and  two  Al(Oh)  near- 
est-neighbors.  The  other  0  is  fourfold-coordinated  in  the  bulk  but 
is  missing  its  Al(Td)  nearest-neighbor  at  the  surface.  The  Al(Td)  is  a 
strong  Lewis-acid  site  on  the  y-Al203  surface  [22,34],  and  there  are 
no  unsaturated  Al(Oh)  sites  on  the  (1 1  l)a  surface.  Upon  relaxing 
the  slab,  with  the  lattice  constants  fixed  but  all  atoms  uncon¬ 
strained,  it  was  found  that  the  atomic  displacements  were  small 
(<0.1  A)  beyond  the  outermost  few  planes,  in  agreement  with  pre¬ 
vious  results  [18]. 

2.3.  Cluster  models  for  y-Al203 

Three  different  cluster  models  (Fig.  2)  were  used,  all  of  which 
are  described  in  detail  elsewhere  [3,4].  All  are  derived  from  the 
(lll)a  defective-spinel  surface  discussed  above.  One  is  an 
AI20O30  structure  with  unsaturated  edge  atoms  [3],  and  the  others 
are  AI8O12  and  AI32O48  clusters  with  pseudo-hydrogen  termination 
(PHT)  [4,35].  A  pseudo-hydrogen  (PH)  is  a  neutral  atom  with  a  non¬ 
integer  nuclear  charge  and  electron  occupancy.  Each  PH  supplies  to 
a  coordinatively-unsaturated  Al  or  0  the  amount  of  charge  that 
would  be  derived  from  the  missing  nearest-neighbor  in  the  y- 
AI2O3  lattice.  Due  to  the  presence  of  both  four  and  sixfold-coordi¬ 
nated  Al  atoms,  four  different  PHs  are  needed  which  are  shown 
in  Fig.  2b  as  H1/2,  etc.  An  H1/2,  for  example,  has  a  nuclear  charge 
of  +1/2  |e|  and  an  electron  occupancy  of  -1/2  |e|.  H1/2  and  H3/4 
are  used  to  terminate  0  atoms,  and  H5/4  and  H3/2  are  used  to 
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Fig.  2.  Cluster  models  used  in  this  study,  (a)  AI20O30  with  unsaturated  edge  atoms, 
(b)  AI8O12  with  PHT,  (c)  AI32O48  with  PHT.  For  clarity  the  AI32O48  pseudo-Hs  are  not 
shown,  and  different  scales  have  been  used  for  the  various  structures.  H1/2,  etc.  are 
pseudo-Hs. 

terminate  A1  atoms.  For  clarity  the  PHs  in  the  AI32O48  model  have 
been  omitted  from  Fig.  2c. 

The  adsorption  energy,  corrected  for  basis  set  superposition  er¬ 
ror  (BSSE),  is  defined  as 

AFads  =  ^(cluster  +  mol)  -  F(cluster)  -  F(mol)  +  AF(BSSE) 

where  the  Es  are,  respectively,  the  relaxed  total  energies  of  the  clus¬ 
ter  with  the  adsorbed  molecule,  of  the  bare  cluster  and  of  the  free 
molecule.  The  last  term  is  the  counterpoise  correction  for  BSSE.  A 
negative  AEads  indicates  an  exothermic  adsorption  process.  Calcula¬ 
tions  for  the  clusters  with  PFIT  were  done  using  ADF  (which  permits 
the  use  of  PHs)  and  the  TZP  basis  sets  described  above.  Eor  technical 
reasons  these  were  done  using  the  (VWN  +  PBE),  rather  than  the 
B3LYP,  functional.  Eor  continuity  with  previous  work  [3],  calcula¬ 


tions  for  the  AI20O30  cluster  (without  PHT)  were  done  using  Gauss¬ 
ian  03  and  the  B3LYP  functional.  This  model  may  be  viewed  as  an 
AI8O12  sub-cluster  situated  within  a  larger  cluster.  The  molecule 
and  the  Al(-0-)3  adsorption  site  were  modeled  with  6-31  lG(d)  ba¬ 
sis  sets,  while  6-31  G(d)  basis  sets  were  used  for  the  rest  of  the 
AI8O12  sub-cluster.  Atoms  outside  the  sub-cluster  were  modeled 
with  3-21 G  basis  sets.  The  AI20O30  AEads  results  obtained  here  for 
Sarin  and  DMMP  (see  below)  differed  by  <1  kcal/mol  from  previous 
results  [3]  obtained  using  6-311G(df)  basis  sets  for  the  molecule 
and  the  Al(-0-)3  site.  The  AI20O30  results  showed  both  internal 
and  spin-restricted/spin-unrestricted  instabilities  [13]  due  to  the 
unsaturated  edge  atoms.  Reoptimizing  the  wavefunction  to  elimi¬ 
nate  these  instabilities  lowered  the  total  energy  (with  or  without 
adsorbates)  by  about  3  kcal/mol.  In  computing  AEads  this  effect 
approximately  cancelled,  leading  to  a  net  reduction  in  adsorption 
energy  of  only  ~0.4  kcal/mol  relative  to  the  uncorrected  result. 

Geometry  optimization  began,  for  the  bare  clusters,  with  relax¬ 
ing  the  positions  of  the  Al(-0-)3  adsorption  site  and  of  the  PHs  (if 
present)  while  keeping  all  other  A1  and  0  atoms  fixed  in  the  ideal- 
lattice  positions.  For  adsorption  the  PHs  remained  fixed  in  the  opti¬ 
mized  positions  while  atoms  in  the  molecule  and  in  the  Al(-0-)3 
adsorption  site  were  allowed  to  vary.  Optimizations  involving 
AI8O12  were  done  using  the  same  convergence  criteria  as  for  the 
free  molecules  (see  above).  For  AI32O48  the  convergence  criteria 
were  relaxed  to  the  ADF  defaults.  For  consistency  in  obtaining  AEads 
the  same  criteria  were  used  for  all  three  components  (adsor¬ 
bate  +  cluster,  bare  cluster  and  free  molecule). 

Vibrational  normal  modes  were  also  obtained  for  adsorbates  on 
the  AI20O30  and  AI8O12  (PHT)  clusters.  Due  to  the  large  number  of 
atoms  involved  this  was  not  attempted  for  the  AI32O48  (PHT)  clus¬ 
ter.  The  interest  here  is  in  the  internal  modes  of  the  molecule,  par¬ 
ticularly  in  v(P=0),  the  frequency  of  the  P=0  stretching  vibration. 
These  receive  little  or  no  contribution  from  displacements  of  atoms 
in  the  cluster  itself  (including  the  PHs).  Hence,  the  results  are  con¬ 
sidered  to  be  reliable  even  though  only  partial  geometry  optimiza¬ 
tions  are  performed,  as  noted  above. 

3.  Results 

3.1.  Free-molecule  structures,  vibrational  frequencies  and  ionization 
potentials 

Table  1  shows  computed  structural  parameters  for  the  free  mol¬ 
ecules  together  with  data  for  TCPO  from  microwave  (p-wave)  spec¬ 
troscopy  [36]  and  for  TMPO  from  X-ray  crystallography  [37].  To  our 
knowledge,  no  corresponding  structural  parameters  are  available 


Table  1 

Observed  and  calculated  structural  and  vibrational  parameters  for  free  molecules^ 


r(F^O) 

r(P-X) 

Z(0=P-X) 

Z(X-P-X) 

v(F^O)'’ 

TCPO'^ 

1.455 

1.989 

d 

103.7 

1321 

1.462  (1.474) 

2.042  (2.030) 

114.9  (114.6) 

103.5  (103.9) 

1268  (1274) 

Sarin 

1311 

1.473  (1.481) 

1256  (1260) 

DMMP 

1276 

1.482  (1.490) 

1223  (1230) 

TMPO^ 

1.489 

1.771 

112.8 

105.9 

1148 

1.498  (1.503) 

1.830  (1.829) 

114.1  (113.9) 

104.4  (104.7) 

1178  (1193) 

^  Observed  values  are  italicized.  Non-italicized  values  are  computed.  The  first  value  is  at  the  B3LYP/6-311G(d)  level,  and  the  value  in  parentheses  is  at  the  (VWN  +  PBE/TZP) 
level.  X  =  Cl  for  TCPO  and  CH3  for  TMPO.  Bond  lengths  are  in  Angstroms,  and  bond  angles  are  in  degrees.  No  experimental  structural  parameters  are  available  for  DMMP  or 
Sarin. 

^  v(P=0)  is  in  cm“\  The  experimental  values  for  TCPO  [43],  DMMP  [44]  and  Sarin  are  from  gas  phase  infrared  spectra.  The  TMPO  value  is  for  the  solid  [45].  The  Sarin  value 
is  an  estimate  [3]  based  on  data  for  the  liquid  [2].  Calculated  values  all  pertain  to  the  gas  phase. 

Experimental  structural  parameters  are  from  Ref.  [36]. 

^  No  experimental  value  is  given. 

®  Experimental  structural  parameters  are  given  in  Ref.  [37].  In  the  crystal  one  C  atom  lies  in  a  mirror  plane  together  with  the  P=0  group  and  is,  therefore,  inequivalent  to 
the  other  C  atoms  which  lie  above  or  below  the  plane.  This  slight  inequivalence  is  ignored  in  the  table,  and  the  experimental  values  given  are  averages  over  all  CH3  groups. 
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Table  2 

Calculated  HOMO  energies  (an)  ionization  potentials  (eV)  for  free  molecules^ 


TCPO 

Sarin 

DMMP 

TMPO 

HOMO  energy  (eh) 

9.51  (8.04) 

8.30  (6.98) 

7.82  (6.53) 

7.19  (5.93) 

Calculated  (vertical  IP) 

11.77  (11.06) 

10.60  (10.17) 

10.17  (9.80) 

9.62  (9.53) 

Calculated  (adiabatic  IP) 

11.25  (10.82) 

9.90  (9.60) 

9.47  (9.29) 

9.10  (9.12) 

Experimental  (vertical  IP) 

11. 89'^ 

11.14° 

10.71*’ 

9.88^ 

^  Values  not  in  parentheses  are  calculated  at  the  B3LYP/6-311G(d)  level.  Values  in  parentheses  are  calculated  at  the  (VWN  +  PBE)/TZP  level. 
From  PES  data,  Ref.  [46]. 

No  experimental  value  is  available.  For  Sarin,  a  VIP  of  11.14  eV  is  estimated  (see  text)  based  on  the  DMMP  results. 

^  From  PES  data.  Ref.  [47]. 


for  either  DMMP  or  Sarin;  although,  moments  of  inertia  are  avail¬ 
able  from  p-wave  data  [38,39].  Generally  good  agreement  is  seen 
between  the  two  different  calculations  and  between  observed 
and  calculated  values  for  TCPO  and  TMPO. 

Sarin  and  DMMP  exhibit  different  conformers  which  have  been 
studied  using  p-wave  spectroscopy  and  ab  initio  or  semi-empirical 
theory  [38-42].  In  the  present  work  the  free  molecules,  upon 
geometry  optimization,  relaxed  into  the  lowest-energy  conforma¬ 
tions  as  reported  previously  [38-40]  and  shown  in  Fig.  1.  For  Sarin, 
calculations  of  the  relaxed  total  energy  vs.  the  angle  of  rotation  of 
the  (CH3)2(H)C-  group  about  the  C-OP  bond  were  done  in  10°  steps 
at  both  the  (WWN  +  PBE)/TZP  and  B3LYP/6-311G(d)  levels.  The  re¬ 
sults  (not  shown)  were  very  similar  to  those  obtained  previously 
[40]  at  the  B3LYP/6-31G(d,p)  level. 

Table  1  also  shows  observed  and  calculated  values  for  v(P=0), 
which  is  sensitive  to  adsorption  (see  below).  The  experimental  val¬ 
ues  for  TCPO  [43],  DMMP  [44]  and  Sarin  (an  estimate  [3]  based  on 
data  for  the  liquid  [2])  are  all  for  the  gas  phase  and  can  be  com¬ 
pared  directly  with  the  computed  results.  On  the  other  hand  TMPO 
is  a  solid,  and  infrared  data  were  obtained  for  a  powder  dispersed 
in  mineral  oil  [45].  v(P=0)  is  known  to  be  significantly  red-shifted 
in  the  condensed  phase  due  to  intermolecular  interactions,  which 
are  not  included  in  the  calculation.  For  example,  v(P=0)  is 
31  cm“^  lower  in  liquid  vs.  gas-phase  DMMP  [44].  Except  for  solid 
TMPO,  the  experimental  results  are  consistently  51  ±  4  cm“^  high¬ 
er  than  the  calculated  values.  This  suggests  a  systematic  error  in 
the  calculation  which  should  cancel,  approximately,  in  obtaining 
Av(P=0),  the  shift  caused  by  adsorption.  It  was  also  found  that 
v(P=0)  for  free  Sarin  is  essentially  independent  of  conformation. 
At  the  B3LYP/6-311'"'"G(d,p)  level  a  difference  of  only  about 
4  cm“^  was  found  between  the  lowest-energy  (“Sarin-l”)  and  high- 
est-energy  (“Sarin-Ill”)  conformations  as  defined  in  Ref.  [40]. 

The  orbital  energies  of  the  molecular  HOMOs  are  needed  in  the 
following  discussion.  Since  these  bear  some  relation  to  experimen- 
tally-observable  IPs  these  quantities  have  been  computed  for  the 
free  molecules.  Table  2  shows  the  VIP  and  AlP  results,  together 
with  the  HOMO  orbital  energies  for  the  neutral  species.  The  VIPs 
can  be  compared  directly  with  PES  data  [46,47]  which  show  better 
agreement  in  the  case  of  B3LYP/6-311G(d).  Eor  B3LYP  the  largest 
discrepancy  between  experiment  and  calculation  (0.54  eV)  is  for 
DMMP,  which  is  chemically  similar  to  Sarin.  Correcting  the  com¬ 
puted  VIP  for  Sarin  (10.60  eV)  by  this  amount  leads  to  an  estimate 
of  about  11.14  eV  for  the  actual  VIP. 

Eor  TCPO  and  TMPO  the  HOMO  is  twofold  degenerate,  with 
large  contributions  from  p;^  and  py  NBOs  on  the  0  atom  of  the 
P=0  group  (which  is  aligned  along  the  z-axis).  Eor  Sarin  and  DMMP 
the  degeneracy  is  lifted  with  splittings  of  about  0.05  and  0.09  eV, 
respectively,  between  the  HOMO  and  HOMO-1.  This  effect  is  ne¬ 
glected  in  Table  2  which  gives  IPs  derived  from  ionization  of  the 
HOMO.  All  quantities  show  the  same  trend,  decreasing  in  magni¬ 
tude  monotonically  in  the  order  TCPO  >  Sarin  >  DMMP  >  TMPO. 
This  reflects  the  inductive  effect  of  strongly-  vs.  weakly-electroneg¬ 
ative  species  (e.g..  Cl  vs.  CH3)  on  the  P  atom  and,  in  turn,  on  the 


HOMO  energy.  It  is  noted  in  passing  that,  in  computing  the  AlP,  a 
Jahn-Teller  distortion  (which  breaks  the  C3V  symmetry)  was  found 
for  the  TCPO  and  TMPO  cations.  Eurther  discussion  of  this  aspect  is 
beyond  the  scope  of  the  present  work. 

3.2.  y-AhOs  Al(Td)  surface  state 

Empty  surface  states  can  be  identified  using  electron  energy 
loss  spectroscopy  (ELS).  As  in  optical  absorption,  the  transition  en¬ 
ergy  includes  the  effect  of  electron-hole  interaction  and  thus  pro¬ 
vides  only  a  lower  limit  on  the  separation  between  the  initial  and 
final  levels  in  the  electronic  ground  state.  Inverse  photoemission 
spectroscopy  [48]  or  tunneling  spectroscopy  can  also  be  used  to 
study  empty  surface  states,  without  the  perturbing  effects  of  core 
or  valence  holes,  but  such  data  are  not  to  our  knowledge  available 
for  y-Al203. 

One  ELS  experiment  was  reported  [49]  for  a  (11  l)-oriented  y- 
AI2O3  film  grown  by  thermal  oxidation  of  a  NiAl  (110)  surface. 
No  surface  states  were  seen  in  the  band  gap.  However,  the  (111) 
surface  in  this  case  is  believed  [50]  to  be  0-terminated  and,  hence, 
would  not  be  expected  to  show  empty  Al-derived  surface  states. 
Ultra-thin  y-Al203  films  have  been  formed  by  oxidation  of  a 
Ni3Al(100)  surface  [51].  These  show  a  band  gap  of  Eg  4.3  eV, 
much  smaller  than  that  for  defect-free  bulk  y-Al203  (Eg  =  7.2  eV 
[52]).  This  was  ascribed  to  the  presence  of  defects,  and  in  any  case, 
such  ultra-thin  oxide  films  may  not  be  representative  of  the  bulk 
material  [53].  Thus  the  experimental  situation  regarding  empty 
surface  states  on  y-Al203  is  uncertain  at  present. 

A  theoretical  study  has  been  performed  [18]  for  various  y-Al203 
surfaces  using  plane-wave  pseudopotentials  and  a  2-DPS  model 
based  on  the  defective-spinel  structure.  The  computed  bulk  band 
gap  (Eg  =  3.97  eV)  underestimates  the  experimental  value,  as  is 
typical  in  LDA  calculations  [54].  On  the  (1 1  l)a  surface  the  valence 
band  maximum  (VBM)  is  found  to  move  down,  and  the  conduction 
band  minimum  (CBM)  to  move  up,  relative  to  the  respective  bulk 
band  edges  which  leads  to  a  surface  band  gap  of  about  4.7  eV.  On 
this  surface,  a  narrow  band  of  states  derived  from  empty  Al(Td) 
orbitals  is  found  at  about  1.3  eV  below  the  surface  CBM.  Similar 
calculations  have  been  done  for  the  (100)  and  (110)  surfaces  of 
the  non-spinel  structure  [22].  These  also  show  distinct  surface- 
state  features  near  the  CBM  which  are  assignable  to  each  of  the  dif¬ 
ferent  unsaturated  A1  sites. 

To  obtain  an  independent  assessment  of  the  surface-state  en¬ 
ergy,  calculations  were  performed  for  the  bulk  lattice  and  for  the 
2-DPS  as  described  above.  Eor  the  bulk  lattice  an  Eg  of  6.2  eV  was 
found  which  is  in  fair  agreement  with  experiment  (Eg  =  7.2  eV 
[52]).  A  previous  B3LYP  treatment  for  an  unspecified  form  of 
AI2O3  (presumably  a-Al203)  gave  a  bulk  Eg  of  8.5  eV  [17]  vs.  the 
experimental  a-Al203  value  of  8.8  eV  [52].  Eig.  3  shows  the  density 
of  states  for  the  2-DPS  in  the  vicinity  of  the  band  edges.  A  surface  Eg 
of  about  7.0  eV  is  found  with  the  lowest-energy  surface-state  (SI) 
located  just  below  (within  ~0.5  eV  from)  the  CBM.  A  second  sur¬ 
face  state  (S2),  also  derived  from  surface  Al(Td)  sites,  is  found  at 
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Fig.  3.  Density  of  states,  in  the  vicinity  of  the  band  edges,  computed  in  the  present 
work  for  the  2-DPS  model  of  Y-AI2O3.  The  bulk  band  gap  and  the  positions  of  the 
empty  surface  states  (SI  and  S2)  are  indicated.  The  zero  of  energy  is  at  the  VBM.  The 
inset  shows  a  schematic  energy  level  diagram  (not  to  scale)  with  definitions  of 
various  quantities. 


about  4.8  eV  above  the  CBM.  The  larger  Eg  for  the  slab  vs.  the  bulk 
lattice  is  consistent  with  previous  results  [18].  The  position  of  SI  is 
essentially  in  agreement  with  that  given  previously  [18];  although, 
the  present  work  finds  it  lying  somewhat  closer  to  the  CBM. 

It  now  remains  to  be  seen  how  closely  the  various  cluster  mod¬ 
els  reproduce  the  2-DPS  SI  position.  The  AI8O12  cluster  (with  PHT) 
gives  a  “band  gap”  of  about  4.4  eV.  This  is  taken  as  the  energy  dif¬ 
ference  between  the  HOMO  and  the  lowest-energy  virtual  state  not 
receiving  any  substantial  contribution  from  atomic  orbitals  local¬ 
ized  on  the  unsaturated  Al(Td)  “surface”  cation.  The  “surface  state” 
associated  with  this  site  then  falls  at  about  1.0  eV  below  the 
“CBM”.  Similar  results  were  found  [35]  for  an  AI8O12  cluster  with 
PHT  designed  to  model  the  ot-Al203  (0001)  surface.  The  PHs  in  this 
case  are  different  from  those  used  for  y-A^Os  since  all  A1  sites  in  a- 
AI2O3  are  Al(Oh),  and  there  are  no  cation  vacancies.  For  this  model 
[35],  after  relaxation.  Eg  ^  4.5  eV  was  found  with  the  surface  state 
at  about  0.5  eV  below  the  “CBM”.  The  surface  state  in  this  case 
arises  from  the  threefold-coordinated  Al(Oh). 

The  AI32O48  cluster  (with  PHT)  gave  Eg  =  4.3  eV  with  the  surface 
state  at  about  0.6  eV  below  the  CBM.  As  shown  in  Fig.  2c,  the  “sur¬ 
face”  of  this  cluster  comprises  seven  Al(Td)  sites.  The  central  Al, 
which  is  the  adsorption  site  in  the  model  calculation,  is  different 
from  the  others  in  that  it  is  threefold-coordinated  to  0  atoms. 
The  other  such  Al  sites  are  all  bonded  to  a  PH  (not  shown  in 
Fig.  2c)  to  compensate  for  a  missing  Al-0  back-bond.  The  SI  state 
is  taken  as  the  virtual  state  in  the  gap  which  is  composed  largely  of 
atomic  orbitals  associated  with  this  central  atom.  Other  such  gap 
states,  derived  from  orbitals  centered  on  the  peripheral  Al(Td)  sites, 
lie  higher  in  energy  and  are  not  involved  in  chemisorption  (see  be¬ 
low).  On  the  (111  )a  surface  of  the  2-DPS  all  Al  sites  are  equivalent 
to  the  central  Al(Td). 

For  the  AI20O30  cluster  (without  PHT)  there  are  many  states  in 
the  “band  gap”  due  to  the  unsaturated  Al  and  0  atoms  at  the 
periphery  of  the  cluster.  The  “bulk  VBM”  was  determined  by  iden¬ 
tifying  the  highest  occupied  molecular  orbital  with  large  contribu¬ 
tions  from  fully-coordinated  0  atoms.  Likewise  the  “bulk  CBM” 
was  identified  by  finding  the  lowest  virtual  orbital  with  large  con¬ 


tributions  from  fully-coordinated  Al  atoms.  The  result  is  Eg  ^ 
5.8  eV  with  the  Al(Td)  surface  state  at  ~0.2  eV  below  the  CBM. 
Hence,  all  three  cluster  models  predict  approximately  the  correct 
location  for  SI  relative  to  the  CBM  in  comparison  to  the  B3LYP  2- 
DPS  result. 

The  question  of  occupied  surface  states  derived  from  0  atoms 
has  also  been  considered.  For  the  (lll)a  surface  none  are  pre¬ 
dicted  in  2-DPS  calculations  to  occur  above  the  bulk  VBM,  either 
here  or  in  previous  work  [18],  and  none  are  observed  experimen¬ 
tally  in  ultraviolet  photoemission  spectroscopy  data  [50,55].  For 
the  AI32O48  PHT  cluster,  however,  states  within  eV  of  the  VBM 
are  derived  largely  from  surface  0  atoms  (i.e.,  those  in  the  upper¬ 
most  0  layer).  This  affects  the  estimates  of  Eg  given  above  but 
not  the  central  issue,  which  is  the  position  of  SI  relative  to  the 
CBM. 

3.3.  Energy  level  alignment 

It  is  now  necessary  to  define  a  common  reference  for  the  energy 
levels  of  the  cluster  and  the  molecule.  Formation  of  the  dative  bond 
(see  below)  involves  the  transfer  of  only  a  small  amount  of  charge 
(~0.1-0.2  |e|)  from  the  P=0  group  to  SI.  Hence,  it  is  meaningful  to 
consider  the  chemical  potential  of  an  electron  in  the  HOMO  or  in 
SI.  Janak  [56]  has  shown  that  Si  =  dEjdni  where  Si  is  the  DFT  energy 
of  orbital  (i)  and  n/  is  the  orbital  population  (treated  as  a  continu¬ 
ous  variable  in  DFT).  For  a  non-metallic  system,  in  which  many- 
body  screening  effects  can  be  neglected,  E  is  effectively  the  total 
energy.  Therefore  in  the  following  discussion  8h,  the  DFT  orbital  en¬ 
ergy  of  the  free-molecule  HOMO,  is  taken  as  the  chemical  potential 
of  an  electron  in  the  HOMO.  Likewise  £5.  the  DFT  orbital  energy  of 
SI,  is  taken  as  the  chemical  potential  of  an  electron  in  the  surface 
state.  The  resulting  £h  and  £5  values  are  summarized  in  Table  3  and 
shown  schematically  in  Fig.  3. 

An  electron  affinity  of 1  eV,  which  places  the  CBM  relative  to 
vacuum  (cf.  Fig.  3),  has  been  estimated  [57]  for  a-Al203.  A  7  of 
1.35  eV  has  been  estimated  [58]  for  y-Al203  on  the  basis  of  electron 
tunneling  measurements  [59]  on  an  ultra-thin  film  on  NiAl(l  1 0) 
[50].  As  noted  above,  such  films  may  not  be  characteristic  of  the 
bulk  material.  If  one  assumes,  nevertheless,  that  x  is  slightly  larger 
than  1  eV  for  y-Al203  and  places  SI  within  ~0.5  eV  of  the  CBM,  as 
found  in  the  present  2-DPS  results,  then  £5  ^  1.5  eV  is  obtained  as 
an  estimate  for  the  actual  value.  This  is  close  to  the  AI32O48  result 
of  1.40  eV.  Hence  one  expects,  a  priori,  that  the  AI32O48  cluster  will 
provide  the  most  accurate  model  of  the  active  site  on  the  defective- 
spinel  (1 1 1  )a  surface,  a  view  which  will  be  supported  by  results  gi¬ 
ven  below. 


Table  3 

Calculated  Afads  for  different  molecules  and  7-AI2O3  model  clusters^ 


^  All  energies  are  kcal/mol  and  have  been  corrected  for  BSSE.  BSSEs  computed  for 
AI8O12  was  also  used  for  AI32O48  since  the  adsorption  geometry  and  computational 
conditions  are  essentially  the  same  in  both  cases. 

^  Results  for  AI20O30  (no  PH  termination)  obtained  at  the  B3LYP/6-311G(d)  level 
and  corrected  for  the  effects  of  wavefunction  instabilities  (see  text). 

HOMO  orbital  energies  calculated  at  the  B3LYP/6-311G(d)  level,  from  Table  2. 

^  Results  for  AI8O12  and  AI32O48  (both  PH-terminated)  obtained  at  the 
(VWN  +  PBE)/TZP  level. 

^  HOMO  orbital  energies  calculated  at  the  (VWN  +  PBE)/TZP  level,  from  Table  2. 

^  Surface  state  (SI )  orbital  energy  determined  as  described  in  the  text.  A  positive 
value  indicates  that  SI  lies  below  the  vacuum  level. 
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3.4.  Adsorption  energies,  geometries  and  vibrational  frequencies 

The  main  theme  of  this  work  is  the  dependence  of  AEads  on  su 
and  as.  The  essential  results  are  summarized  in  Table  3.  Fig.  4a 
shows,  as  an  example,  the  optimized  structure  of  TMPO  adsorbed 
on  the  Al80i2  cluster  (with  the  PHs  removed  for  clarity).  For  other 
adsorbates  the  structure  is  similar  in  appearance  (e.g..  Ref.  [3]).  For 
a  given  cluster,  comparing  results  for  different  molecules  shows 
the  correlation  between  AFads  and  su.  As  one  expects  intuitively 
on  the  basis  of  the  frontier-orbital  concept  [7-9],  decreasing  an 
makes  charge  transfer  easier  and  leads  to  an  increased  AFads-  For 
a  given  molecule,  comparing  results  for  the  AI8O12  vs.  AI32O48  clus¬ 
ters  shows  the  correlation  between  AFads  and  as.  Again  as  expected, 
moving  SI  up  toward  vacuum  makes  the  surface  less  acidic  (i.e.,  a 
weaker  acceptor)  and  impedes  dative-bond  formation,  thus 
decreasing  the  magnitude  of  AFads- 

Fig.  5  shows  plots  of  the  results  in  Table  3.  All  clusters  yield  an 
essentially  linear  dependence  of  AFads  on  (an-as)  with  slightly  dif¬ 
ferent  parameters  for  AI20O30  vs.  AI8O12  and  AI32O48.  Initially  the 
AI8O12  and  AI32O48  results  were  fitted  separately,  but  the  fits  were 
virtually  indistinguishable  and  Fig.  5  shows  a  single  fit  for  both. 
The  reason  for  the  small  differences  between  the  PHT  and  non- 
PHT  parameters  (about  5%  in  the  slope  and  2%  in  the  intercept)  is 
difficult  to  identify.  One  possibility  is  the  different  DFT  functionals 
used  for  the  two  sets  of  calculations. 

Although  the  qualitative  dependence  of  AFads  on  an  and  as  is  ex¬ 
pected,  the  quantitatively  linear  behavior  in  Fig.  5  is  not.  A  linear 
combination  of  atomic  orbitals  (LCAO)  treatment  (e.g..  Refs.  [7,9]) 
yields  AFads  =  H^/Aa  where  H  is  the  interaction  energy  between 
the  empty  SI  state  and  the  doubly-occupied  molecular  HOMO 
and  Aa  =  as-aH.  This  is  based  on  a  perturbation  approximation 
(|H|  <  |Aa|)  which  does  not  apply  in  the  present  case  since,  for 
the  AFads  and  Aa  values  encountered  here,  H  and  Aa  would  be  of 


Fig.  4.  Model  showing  the  optimized  structure  for  (a)  TMPO  and  (b)  H2O  adsorbed 
on  the  AI8O12  cluster.  The  pseudo-Hs  have  been  removed  for  clarity.  In  (b)  the 
dashed  line  shows  a  hydrogen  bond  with  r(O-H)  =  1.964  A. 


Fig.  5.  Graphical  display  of  the  results  in  Table  3.  The  linear  regressions  are  given  by 
AFads  = -106.2  +  12.22(£h-£s)  for  AI8O12  (circles)  and  AI32O48  (squares),  r=  0.984 
AFads  =  -104.1  +  12.82(8H-es)  for  AI20O30  (triangles),  r  =  0.995  with  AFads  in  kcal/mol 
and  (sh-ss)  in  eV.  The  r  values  are  the  correlation  coefficients.  Data  points  are 
labeled  for  AI8O12  and  appear  in  the  same  order  for  the  other  plots.  The  open 
symbols  show  AFads  values  computed  for  CO  (see  text)  which  were  not  included  in 
the  linear  regression  analysis. 

similar  magnitude.  Fits  to  AFads  which  were  linear  in  l/As  (not 
shown)  were  not  quite  as  good  as  those  in  Fig.  5.  Different  fits  were 
obtained  for  AI8O12  and  AI32O48  when  plotting  vs.  l/As,  with  corre¬ 
lation  coefficients  of  r  =  0.931,  0.967  and  0.981  for  AI8O12,  AI20O30 
and  AI32O48,  respectively. 

When  the  effects  of  su  and  £5  are  “factored  out”,  similar  behavior 
is  seen  for  different  cluster  models,  functionals  and  basis  sets.  With 
all  else  equal,  differences  in  AFads  for  various  clusters  can  be  traced 
primarily  to  differences  in  £5.  Previous  2-DPS  results  for  y-A103 
[22],  obtained  for  the  non-spinel  bulk  structure,  have  shown  that 
the  Lewis  acidities  of  the  various  types  of  A1  surface  sites  can  be 
correlated  directly  with  the  energies  of  the  corresponding  empty 
surface  states  relative  to  vacuum  (£5  in  the  present  notation).  Ener¬ 
gies  in  the  range  of  -2.5  ^  £s  ^  +0.1  eV  were  reported  (using  a  sign 
convention  opposite  to  that  used  here).  In  particular,  £5  =  -2.5  eV 
was  found  for  the  threefold-coordinated  Al(Td)  site  on  the  (110) 
surface  which  is  close  to  the  value  of  (-)2.18  eV  obtained  here  (Ta¬ 
ble  3)  for  the  AI8O12  (PHT)  cluster.  It  was  also  found  that  the 
strength  of  the  interaction  of  CO  with  these  sites,  as  measured  by 
AFads  and  by  Av(C=0),  correlates  with  the  difference  between  £5 
and  £h  (the  energy  of  the  CO  HOMO,  i.e.,  the  5a  orbital).  A  further 
discussion  of  these  points,  including  CO  adsorption,  is  given  below. 

All  adsorbates  except  TCPO  show  the  possibility  of  H-bonding 
between  H  atoms  in  the  alkyl  group(s)  and  0  atoms  in  the  AI2O3 
cluster.  This  is  based  on  the  observation,  in  the  optimized  struc¬ 
tures,  of  one  or  more  CH-0  distances  in  the  range  of  2.1 -2.4  A. 
Similar  effects  have  been  noted  in  other  studies  [5,6]  of  CWA 
adsorption  on  ionic  oxides.  No  attempt  has  been  made  here  to  as¬ 
sess  the  contribution  of  these  effects  to  AFads-  However,  ab  initio 
studies  of  other  systems  [60,61]  suggest  that  such  CH-0  bonds 
are  weak,  with  a  bond  energy  of  at  most  ~3  kcal/mol. 

Cluster  models  of  the  sort  used  here  neglect  the  long-range  con¬ 
tribution  to  the  electrostatic  potential  at  the  adsorption  site  that 
would  arise  in  a  semi-infinite  lattice.  The  possibility  was  consid¬ 
ered  that  a  difference  in  this  potential  could  contribute  to  the  dif- 
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ference  in  AEads  for  the  various  clusters.  With  the  active  Al(Td)  as 
the  test  site  (J\  tho  sum  of  gz/ry  over  all  other  sites  (i)  including 
the  PHs  gave  -1.232  and  -1.173  |e|/A  for  AI8O12  and  AI32O48, 
respectively.  Here  g/  is  the  ionic  charge  computed  using  the  multi¬ 
pole-derived  charge  (quadrupole),  or  MDC-Q,  method  of  Swart 
et  al.  [62],  and  Vy  is  the  distance  from  (j).  Likewise  the  MDC-Q. 
charge  on  the  active  Al(Td)  site  of  the  bare  cluster  is  similar  for 
AI8O12  (+1.82)  vs.  AI32O48  (+1.77).  For  comparison,  a  Mulliken 
charge  of  +1.66  is  found  for  the  surface  Al(Td)  ion  in  the  2-DPS  cal¬ 
culation  described  above.  For  the  AI20O30  cluster,  a  Mulliken  charge 
of  +1.30  is  found  for  the  “surface”  Al(Td),  which  is  smaller  than  the 
corresponding  values  given  above.  However,  the  sum  over  all  qilry 
is  -1.231  |e|/A,  essentially  the  same  as  that  found  for  the  other 
clusters.  Hence,  the  electrostatic  environment  in  the  vicinity  of 
the  adsorption  site  is  similar  for  all  three  clusters;  although,  it 
may  of  course  differ  from  that  for  a  semi-infinite  crystal. 

Examination  of  the  cluster  orbitals  shows  that  adsorption  re¬ 
moves  SI  from  the  gap.  Bonding  orbitals,  with  large  components 
from  the  phosphonyl  0  atom  and  small  contributions  from  the 
Al(Td)  adsorption  site,  appear  a  few  eV  below  the  “VBM”.  This  is  di¬ 
rect  (theoretical)  evidence  for  the  involvement  of  SI  in  the  chemi¬ 
sorption  process.  The  MDC-Q  results  for  AI8O12  and  AI32O48  show 
that  the  charge  transferred  to  the  Al  atom  is  the  range  of  -0.11 
to  -0.21  |e|,  depending  on  the  cluster  and  the  adsorbate. 

The  configuration  of  adsorbed  Sarin  differs  slightly  from  the 
lowest-energy  Sarin-1  structure  [40]  of  the  gas-phase  molecule. 
For  the  AI20O30  cluster  the  isopropyl  group  rotates  about  the  C- 
OP  bond  by  about  27°  in  the  direction  of  Sarin-11;  whereas,  for 
the  PHT  clusters  the  rotation  is  by  about  13°  in  the  direction  of  Sar- 
in-111.  For  all  models.  Table  4  shows  that  r(Al-OP)  decreases  with 
increasing  AEads  as  expected.  All  models  yield  about  the  same 
(A1-0=P)  angle  for  a  given  molecule  except  for  Sarin  where  a  lar¬ 
ger  angle  is  found  for  AI20O30  (153°)  vs.  AI8O12  and  AI32O48  (131°). 
This  probably  results  from  the  different  Sarin  configurations,  noted 
above,  for  AI20O30  vs.  AI8O12  and  AI32O48.  The  orientation  of  the  iso¬ 
propyl  group  for  Sarin  adsorbed  on  AI20O30  requires  a  larger  (Al- 
0=P)  angle  to  avoid  steric  interaction  between  one  of  the  CH3 
groups  and  the  cluster. 

The  computed  Av(P=0)  values  for  AI20O30  and  AI8O12  are  in 
fairly  good  agreement  with  each  other  and  (where  available)  with 
experimental  data  except,  again,  in  the  case  of  Sarin  for  which  the 
AI8O12  result  is  too  large  vs.  experiment  (-66cm“^  [2,3]).  The 
explanation  for  this  is  uncertain  at  present.  It  may  lie  in  the  differ¬ 
ent  Sarin  conformations  in  the  two  models;  however,  as  noted 
above,  v(P=0)  for  free  Sarin  appears  to  be  essentially  independent 


of  conformation.  For  adsorbed  DMMP  the  computed  P=0  stretch  is 
strongly  coupled  to  modes  at  1161  and  1174  cm“\  and  v(P=0)  is 
taken  as  the  average  of  these  values. 

Results  for  the  gas-phase  molecules  (Table  1)  show  the  ex¬ 
pected  decrease  in  v(P=0)  with  decreasing  bond  strength  (i.e.,  with 
increasing  r(P=0)).  Table  4  shows  that  r(P=0)  is  in  all  cases  greater 
for  the  adsorbed  vs.  the  gas  phase,  consistent  with  a  weakening  of 
the  P=0  bond.  The  increase  in  bond  length  is  small,  being  in  the 
range  of  5r(P=0)  ?^0.03  to  0.06  A,  but  there  does  appear  to  be  a 
rough  correlation  between  Av(P=0)  and  5r(P=0).  However,  the 
correlation  between  Av(P=0)  and  AEads  is  somewhat  more  tenu¬ 
ous,  as  is  seen  by  comparing  results  for  TCPO  and  TMPO.  These 
show  essentially  identical  Av(P=0)  values  but  AEads  values  differ¬ 
ing  by  nearly  a  factor  of  two.  Similar  behavior  was  found  previ¬ 
ously  [63]  for  the  adsorption  of  TCPO  and  DMMP  on  amorphous 
silica  via  H-bonding  between  Si-OH  groups  and  the  phosphonyl 
0  atom.  A  AEads  of  -10.4  (-20.0)  kcal/mol  was  computed  for  TCPO 
(DMMP),  but  the  experimental  Av(P=0)  values  are  about  -34  and 
-19  cm“\  respectively.  The  reason  for  this  behavior  lies  in  a  un¬ 
ique  property  of  TCPO,  whereby  Cl  orbitals  are  strongly  mixed  with 
both  the  HOMO  and  the  P=0  n  orbital.  For  the  other  species  stud¬ 
ied  here  the  HOMO  consists  almost  entirely  of  orbitals  on  the  0 
atom  of  the  P=0  group.  For  TCPO,  the  HOMO  and  the  P=0  ji-bond 
are  in  effect  coupled,  and  perturbation  of  the  former  by  adsorption 
then  leads  to  a  disproportionately  strong  effect  on  the  latter  and  a 
correspondingly  large  Av(P=0).  Hence,  some  caution  must  be  ap¬ 
plied  when  interpreting  Av(P=0)  strictly  in  terms  of  the  strength 
of  the  adsorption  bond. 

Finally,  other  possible  modes  of  adsorption  were  considered, 
even  though  experimental  data  [1,2]  for  DMMP  and  Sarin  clearly 
indicate  A1-0=P  dative  bond  formation  as  the  energetically-fa¬ 
vored  process.  Adsorption  of  Sarin  on  the  AI8O12  cluster  gave  AEads 
of  -26.2  or  -31.0  kcal/mol,  respectively,  for  bonding  between  Al 
and  the  F  or  the  isopropoxy  0  atom.  These  have  not  been  corrected 
for  BSSE,  and  doing  so  would  reduce  the  magnitudes  by  ~2  kcal/ 
mol.  These  can  be  compared  with  the  BSSE-corrected  value  of 
-51.5  kcal/mol  (Table  3)  for  A1-0=P  bonding.  These  results  are 
consistent  with  previous  results  [3]  for  the  AI20O30  cluster.  Bonding 
of  TCPO  to  AI8O12  via  an  Al-Cl-P  bond  was  also  tested.  During 
relaxation  the  molecule  moved  away  from  the  surface,  and  no 
bond  formed.  Eor  TCPO  a  “reverse”  dative  bond  was  investigated, 
wherein  an  0  atom  on  the  AI2O3  surface  acts  as  an  electron  donor 
and  the  P  atom  acts  as  an  acceptor.  This  might,  in  principle,  be 
energetically  feasible  since  the  lowest  unoccupied  molecular  orbi¬ 
tal  of  gas-phase  TCPO  is  computed  to  lie  at  3.485  eV  below  vacuum 


Table  4 

Results  for  adsorption  on  7-AI2O3  model  clusters^ 


^  In  all  cases,  adsorption  is  via  A1-0=P  dative  bond  formation  (cf.  Fig.  4).  Energies  are  in  kcal/mol,  bond  lengths  in  Angstroms,  bond  angles  in  degrees  and  frequencies  in 
cm“L  The  AEads  values  repeat  those  in  Table  3. 

^  Numbers  in  parentheses  are  P=0  bond  lengths  computed  for  the  free  molecules. 

°  Av(P=0)  is  the  calculated  v(P=0)  for  the  adsorbed  species  minus  the  calculated  gas-phase  value.  No  experimental  Av(P=0)  data  are  available  for  TCPO  or  TMPO.  For 
DMMP  the  experimental  result  is  -60  cm“^  (Ref.  [1]).  For  Sarin  the  estimated  experimental  result  is  -66  cm“^  (Refs.  [2,3]). 
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at  the  (VWN  +  PBE)/TZP  level.  The  molecule  was  initially  posi¬ 
tioned  with  P=0  bond  pointing  away  from  the  surface  and  with 
the  three  Cl  atoms  coplanar  with  the  P  atom.  The  AI32O48  cluster 
was  used,  and  the  P  was  positioned  above  a  coordinatively-unsat- 
urated  0  atom  (see  above).  However,  here  again  the  molecule  re¬ 
laxed  away  from  the  surface  with  no  bond  formation. 

3.5.  Comparison  of  2-DPS  and  cluster  results 

There  are,  to  our  knowledge,  no  experimental  data  giving  ATads 
for  the  systems  of  interest  here.  It  is  therefore  difficult  to  assess  the 
absolute  accuracy  of  the  results.  However,  comparison  can  be 
made  with  results  for  a  2-DPS  which  is  assumed  to  model  the  sur¬ 
face  of  a  semi-infinite  crystal  more  accurately  than  does  a  finite 
cluster.  One  such  comparison  is  for  the  physisorption  of  molecular 
H2O  using  the  same  functionals,  basis  sets,  etc.  as  for  the  phospho- 
nyl  species.  In  all  cases  a  dative  bond  (H2O-AI)  forms  between  the 
0  atom  and  the  unsaturated  Al(Td)  site.  A  2-DPS  result  of 
-28.1  kcal/mol  for  the  (1  1  l)a  surface  has  been  found  by  Pinto 
and  Elliott  [64]  using  a  plane-wave  pseudopotential  approach. 
The  results  obtained  here  (corrected  for  BSSE)  are  AEads  =  -40.1, 
-36.5  and  -30.8  kcal/mol,  respectively,  for  the  AI20O30,  AI8O12 
(PHT)  and  AI32O48  (PHT)  clusters.  Again,  AEads  decreases  in  magni¬ 
tude  as  8s  shifts  toward  vacuum.  The  AI32O48  result  is  close  to  that 
found  for  the  2-DPS  which  suggests  that  this  cluster,  of  the  three 
used  here,  gives  the  most  reliable  absolute  values  for  AEads  in  com¬ 
parison  to  the  (1 1  l)a  defective-spinel  surface. 

In  the  case  of  H2O  adsorption  there  is  evidence  for  a  contribu¬ 
tion  to  AEads  from  hydrogen  bonding  between  H  and  surface  0 
atoms.  Pig.  4b  shows,  as  an  example,  the  optimized  structure  for 
the  AI8O12  (PHT)  cluster  which  exhibits  an  H-bond  with  r(0- 
H)  =  1.964  A.  This  complicates  somewhat  an  interpretation  of  AEads 
strictly  in  terms  of  dative-bond  formation.  Nevertheless,  the  trend 
in  AEads  vs.  ss  is  in  accord  with  that  expected  on  the  basis  of  the  re¬ 
sults  for  the  phosphonyl  species. 

Another  comparison  can  be  made  for  the  adsorption  of  CO,  a 
weak  Lewis  base  often  used  to  probe  the  acidity  of  oxide  surfaces. 
Adsorption  occurs  (e.g..  Ref.  [22])  via  an  A1-C=0  dative  bond 
involving  the  5a  orbital  of  CO  which  is  the  HOMO.  A  DPT  calcula¬ 
tion,  as  described  above,  was  done  for  the  (1 1  l)a-oriented  2-DPS 
using  6-31  lG(d)  basis  sets^  for  CO  and  a  (1  x  1)  surface  unit  cell 
with  a  CO  at  every  surface  Al(Td)  site.  The  y-Al203  slab  was  modeled 
using  the  same  basis  sets  as  those  described  above.  It  is  assumed  that 
any  interaction  between  adjacent  molecules  can  be  neglected.  Based 
on  results  [65]  for  CO  adsorbed  on  MgO(lOO)  and  on  the  present 
intermolecular  distance  of  5.6  A,  lateral  interaction  is  expected  to  de¬ 
crease  the  magnitude  of  AEads  by  <0.5  kcal/mol  relative  to  the  value 
for  an  isolated  molecule.  CO  was  placed  identically  on  both  the 
upper  and  lower  slab  faces  in  order  to  maintain  symmetry. 

The  2-DPS  calculation  gave  AEads  =  -10.1  kcal/mol  with  opti¬ 
mized  bond  lengths  of  r(Al-C)  =  2.200  A  and  r{C=0)  =  1.122  A.  Por 
the  free  molecule,  r(C=0)  =  1.127  A  was  found  using  the  same 
(B3LYP)  functional  and  basis  sets.  The  5a  orbital  has  a  slight  anti¬ 
bonding  character;  hence,  adsorption  leads  to  a  very  small  decrease 
in  r(C=0).  The  AI32O48  PHT  cluster  gave  a  stronger  interaction,  with 
AEads  =  -15.7  kcal/mol,  r(Al-C)  =  2.144  A  and  r(C=0)  =  1.133  A.  The 
corresponding  free-molecule  r(C=0)  was  found  to  be  1.141  A  using 
the  VWN  +  PBE  functional  and  TZP  basis  sets.  The  AI8O12  PHT  cluster 
gave  a  still-larger  AEads  of  -19.8  kcal/mol  due  to  the  lower-lying  SI 
(Table  3).  Similar  calculations  for  the  AI20O30,  as  described  above. 


^  The  C  and  0  basis  sets  were  obtained  from  the  Extensible  Computational 
Chemistry  Environment  Basis  Set  Database,  Version  02/02/06,  developed  and 
distributed  by  the  Molecular  Science  Computing  Eacility;  Environmental  and 
Molecular  Sciences  Laboratory  which  is  part  of  the  Pacific  Northwest  Laboratory: 
P.O.  Box  999;  Richland,  WA  99352;  USA. 


gave  AEads  = -20.8  kcal/mol  with  r(Al-C)  =  2.136  A  and 
r(C=0)  =  1.1 17  A.  All  energies  have  been  corrected  for  BSSE. 

Por  free  CO,  an  8h  of  9.015  eV  is  found  at  the  (VWN  +  PBE)/TZP 
level  and  10.360  eV  at  the  B3LYP/6-311G(d)  level.  With  these  re¬ 
sults  and  the  ss  values  given  in  Table  3,  the  calculated  AEads  values 
for  CO  can  be  placed  in  Pig.  5.  In  all  cases  the  results  are  in  good 
accord  with  extrapolations  of  the  linear  fits  to  the  phosphonyl  re¬ 
sults.  As  was  seen  for  H2O,  the  AI32O48  result  gives  the  best  agree¬ 
ment  with,  but  is  larger  than,  the  2-DPS  value  of  AEads  =  -10.1  kcal/ 
mol  for  the  (1 1 1  )a  defective-spinel  surface.  Por  comparison,  2-DPS 
calculations  for  CO  adsorption  at  the  Al(Td)  site  on  the  (110)  sur¬ 
face  of  the  non-spinel  model  give  AEads  =  -18.2  [22]  or 
-17.1  kcal/mol  [33].  This  site  has  a  computed  ss  of  2.5  eV,  as  noted 
above.  It  is  thus  more  acidic  than  the  Al(Td)  site  on  the  AI32O48  clus¬ 
ter,  with  8s  =  1.40  eV,  but  closer  in  acidity  to  that  on  the  AI8O12, 
with  8s  =  2.18  eV.  Thus  the  non-spinel  values  for  AEads  [22,33]  are 
close  to  the  AI8O12  result  of  -19.8  kcal/mol. 


4.  Discussion 

The  present  work  analyzes  a  basic  principle  affecting  dative- 
bond  formation  between  phosphonyl  reagents  and  Lewis  acid  sites 
on  the  y-Al203  surface.  Presumably  the  same  considerations  affect 
similar  adsorption  processes  for  other  systems.  Table  3  and  Pig.  5 
demonstrate  the  sensitivity  of  AEads  to  8h  and  to  8s-  In  the  absence 
of  experimental  data  for  AEads  it  is  necessary  to  compare  the  cluster 
results  with  those  for  a  2-DPS.  Such  comparisons  for  H2O  and  CO 
(see  above)  indicate  that  the  AI32O48  PHT  cluster  at  the 
(VWN  +  PBE)/TZP  level  gives  the  most  accurate  AEads  in  compari¬ 
son  to  the  2-DPS  model  of  the  defective-spinel  (lll)a  surface. 
Therefore  by  implication  this  model,  of  the  three  studied  here, 
gives  the  most  reliable  representation  of  SI  for  that  surface.  On 
the  other  hand,  85  for  the  AI8O12  PHT  cluster  is  close  to  that  found 
for  the  Al(Td)  site  on  the  (110)  non-spinel  surface.  Consistent  with 
this,  the  AEads  values  for  CO  adsorption  are  also  found  to  be  close 
for  these  two  substrates. 

It  might  appear  that  issues  involved  in  the  formulation  of  clus¬ 
ter  models  could  be  avoided  simply  by  working  exclusively  with  2- 
DPS  models.  However,  the  molecules  of  interest  here  and  in  other 
CWA  studies  are  comparable  in  size  to,  or  larger  than,  the  y-Al203 
(1x1)  surface  unit  cell.  Avoiding  significant  steric  interaction  be¬ 
tween  nearest-neighbor  adsorbates  in  a  2-DPS  calculation  requires 
the  use  of  a  suitably  large  slab  supercell  which  may  make  the  cal¬ 
culation  computationally  intractable.  Hence,  clusters  are  still  pref¬ 
erable  (or  even  mandatory),  provided  that  8s  can  be  reliably 
obtained.  Having  constructed  such  a  cluster  model  one  can  then 
proceed  with  some  confidence  that  reasonable  adsorption  results 
will  be  computed. 

It  is  worthwhile  to  re-examine  previous  work  [3]  in  light  of  the 
present  results.  Absolute  adsorption  energies  may  be  difficult  to 
compute  reliably  for  any  model  (cluster  or  2-DPS)  because  of  the 
need  to  reproduce  the  experimental  8s  which  may  itself  be  uncer¬ 
tain,  as  is  the  case  here.  However,  other  quantities  of  interest  in¬ 
clude  the  adsorption  geometry  (i.e.,  lengths  and  angles  of  the 
chemisorption  bonds),  shifts  in  adsorbate  vibrational  frequencies 
relative  to  the  gas  phase  and,  especially,  relative  adsorption  ener¬ 
gies  for  similar  molecules  or  for  different  functional  groups  of 
the  same  molecule.  These  kinds  of  information  can  often  be  ob¬ 
tained  even  in  calculations  that  are  only  qualitative  with  regard 
to  the  absolute  AEads- 

Consider,  for  example,  AEads  for  the  adsorption  of  DMMP  vs. 
Sarin  (Table  3).  The  AI32O48  results  are  in  qualitative  agreement 
with  those  obtained  previously  [3]  for  the  much  simpler  AI20O30 
model  which  show  a  larger  AEads  for  DMMP.  However,  the  differ¬ 
ence  of  1.6  kcal/mol  is  less  than  the  previous  estimate  of 
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8.2  kcal/mol.  The  difference  in  AEads  for  the  two  molecules  is 
now  understood  to  result  from  the  somewhat  smaller  su  of 
DMMP.  Another  example  concerns  adsorption  via  an  A1  bond 
to  the  phosphonyl  vs.  the  alkoxy  0  atom  of  DMMP  or  Sarin. 
For  Sarin  at  the  B3LYP/6-311G(d)  level,  the  simple  AI20O30  clus¬ 
ter  gives  AFads  =  -30.9  kcal/mol  for  bonding  to  the  alkoxy  0 
atom  vs.  -49.0  kcal/mol  for  bonding  to  the  P=0,  a  difference 
of  about  18  kcal/mol.  A  similar  difference,  -20.9  vs.  -39.8  kcal/ 
mol,  is  found  at  the  (VWN  +  PBE)/TZP  level  for  adsorption  on 
the  AI32O48  PHT  cluster. 

When  considering  the  effects  of  different  ligands  on  the  relative 
adsorption  energies  of  phosphonyl  species  it  should  be  possible  to 
use  the  easily-computed  su  of  the  free  molecule  as  a  quantitative 
indicator.  In  view  of  the  above  discussion  concerning  H2O  adsorp¬ 
tion,  such  use  of  8h  involves  an  assumption  that  other  forms  of 
interaction,  such  as  H-bonding,  can  be  neglected.  Likewise,  in  com¬ 
paring  relative  adsorption  energies  for  Lewis-acid  sites  on  other 
oxides,  insight  can  be  obtained  by  computing  (or  measuring)  £5- 
One  assumes  here  that  the  substrate  orbital  involved  in  the  adsorp¬ 
tion  corresponds  to  a  well-defined  surface  state  lying  in  the  band 
gap  and  that  any  modification  of  the  surface  does  not  cause  a  sig¬ 
nificant  change  in  the  electronic  structure  of  the  surface  state. 

Similar  considerations  might  apply  even  when  the  HOMO  is  not 
derived  from  NBOs  on  the  P=0  group.  For  example,  other  work  [4] 
shows  that  adsorption  of  the  nerve  agent  VX  on  OH-free  y-Al203 
also  involves  A1-0=P  bonding.  However,  the  HOMO  and  HOMO- 
1  in  this  case  are  derived  largely  from  NBOs  on  the  nitrogen  and 
sulfur  atoms,  respectively,  neither  of  which  is  directly  involved  in 
chemisorption.  The  orbital  energies  of  the  deeper-lying  NBOs  of 
the  phosphonyl  group  (HOMO-2  and  HOMO-3)  are  computed  to 
be  6.36  and  6.72  eV,  respectively.  Using  the  average  for  £h  in  the 
linear  relationship  obtained  in  Fig.  5  then  gives  AFads  =  -43.4  kcal/ 
mol  for  adsorption  on  AI32O48.  This  compares  reasonably  well  with 
the  calculated  value  [4]  of  -39.2  kcal/mol. 

5.  Conclusions 

Density  functional  theory  has  been  applied  to  a  study  of  the 
adsorption  of  a  series  of  phosphonyl  compounds  on  y-Al203.  The 
results  are  as  follows. 

(1 )  All  the  reagents  studied  adsorb  via  an  A1-0=P  dative  bond  at 
a  coordinatively-unsaturated  Al(Td)  site.  Other  possible 
modes  of  adsorption  are  energetically  less  favorable. 

(2)  The  adsorption  energy  (AFads)  shows  a  linear  dependence 
on  the  difference  between  the  one-electron  orbital  ener¬ 
gies  of  the  molecular  HOMO  (an)  and  the  Al(Td)  surface 
state  (8s).  Trends  in  AFads  for  different  reagents  on  a  given 
cluster,  and  for  the  same  molecule  on  different  clusters, 
can  be  understood  quantitatively  in  terms  of  variations 
in  8h  and  85. 

(3)  Cluster  models  can  be  constructed  which  give  8s  in  good 
agreement  with  two-dimensionally-periodic  slab  (2-DPS) 
predictions  based  on  fundamentally-different  models  for 
the  y-Al203  bulk  lattice.  For  test  molecules  such  as  H2O 
and  CO,  properly-constructed  clusters  then  give  AFads  results 
in  good  agreement  with  those  from  the  corresponding  2-DPS 
calculations.  Cluster  models  can  thus  be  a  viable  alternative 
to  2-DPS  calculations  in  the  treatment  of  adsorption  on  a 
complex  ionic  oxide  such  as  y-Al203. 

(4)  The  results  may  be  extensible  to  other,  non-phosphonyl  spe¬ 
cies.  For  adsorption  of  CO  via  an  A1-C=0  bond,  AFads  values 
obtained  from  2-DPS  calculations  agree  well  with  cluster- 
model  predictions  based  on  a  linear  extrapolation  of  the 
phosphonyl  results.  Thus  AFads  for  phosphonyls,  and  perhaps 


for  other  species  that  adsorb  by  a  similar  dative-bond  forma¬ 
tion,  can  be  estimated  quantitatively  using  the  easily-com¬ 
puted  £h  of  the  free  molecule. 
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